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ABSTRACT 


The high 0 + - H + transition level at latitudes above 50 degrees and the 
plasmapause with its characteristic density plateau suggest that both features 
are strongly correlated; upward fluxes of protons can consistently account for 
both observations. It is shown that plasma escape can produce sufficiently 
large fluxes at latitudes where the field lines are open, thus significantly de- 
creasing the hydrogen-ion density. As a result, plasma diffusion across field 
lines is induced- coupled with upward fluxes of ionization along field lines, which 
leads to the formation of a density plateau that characterizes the transition region 
from the open to the closed magnetosphere. Collision- type perpendicular diffusion 
would produce an enormously rapid transition over a distance of a few meters, 
which is not observed. This suggests that a turbulent state is built up in the 
transition region; thus, Bohm’s coefficient for turbulent diffusion is em- 
ployed, and this leads to a quantitatively satisfying description of the plasma- 
pause. 
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The Plasmapause and its Relation to the Ion 
Composition in the Topside Ionosphere 
H. G. Mayr 

INTRODUCTION 

Whistler data [Carpenter , 1963, 1966] and ion-spectrometer measure- 
ments [Taylor et al , , 1965] have shown that the electron concentration decreases 
abruptly and very rapidly above L values of 3 to 5. Recent observations from 
OGO satellite [Taylor, private communication] confirm this feature of the mag- 
netosphere, which is referred to as the "plasmapause. " 

Mayr, Brace, and Dunham [1967] studied ion composition observed by 
Thomas et al . [1966] for the two major ion constituents 0+ and H + during mini- 
mum solar activity. They found that up to 50 degrees geomagnetic latitude the 
observed latitudinal temperature increase [Brace, Reddy, and Mayr , 1967] 
accounts for the increasing 0 + - H + transition level. At higher latitudes, how- 
ever, the altitude of the observed transition level continues to increase; in con- 
trast to theory, that predicts a slight decrease of this level. In order to explain 
this discrepancy, Mayr, Brace, and Dunham [1967] postulated that there are 
upward fluxes of ionization into the protonospheric field tubes that increase 
with latitude and thus significantly reduce proton concentration. In a static 
model such fluxes can arise only from ionization leakage out of the field tubes 
by diffusion across field lines. It was assumed that this was due to a density 
gradient across the magnetic field, the gradient being maintained through escape 
of protons at higher latitudes where the field line configuration permits it. This 
suggested that the density gradient across high-latitude field lines is related to 
the density plateau that constitutes the plasmapause. A consistent concept was 
thus proposed, combining the plasmapause at high altitudes with ion composi- 
tional features at high latitudes. 
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This paper investigates this concept. Turbulent diffusion postulated by 
Bohm [1949], when combined with plasma escape, is shown to account quantita- 
tively for the plasmapause and its relation to ion composition at high latitudes. 

PHYSICAL CONCEPT 


General 

Describing the magnetosphere and considering the density structure of 
the dominant ion species, proton, we can distinguish three different regions in 
order of increasing altitude. 

(1) The topside ionosphere. During minimum solar activity, the major 
ion constituents in this region (oxygen ions) are almost in diffusive equilibrium. 
The minor constituents (hydrogen ions) are in chemical equilibrium with 0 + , O, 
and H through the resonance charge exchange reaction H + 0 + il! H + + O. 

(2) The transition region between ionosphere and protonosphere. Here, 
diffusion processes become important for protons and comparable with chemical 
processes. In this most complex region, the H + concentration is very sensitive 
to the proton diffusion flux. As shown by Geisler and Bowhill [1965], upward 
fluxes of ionization, if in the order of the critical flux, can decrease the proton 
concentration by almost two orders of magnitude and correspondingly increase 
the transition level. Similarly, downward fluxes can significantly increase the 
H + concentration and lower the transition level. 

(3) The protonosphere. Here, protons are the major constituents and can 
be considered in a state of fast diffusion along field lines, that is almost equiva- 
lent to diffusive equilibrium. The protonosphere has its roots in the previously 
described transition region. Any effect on the hydrogen-ion density there, due 
to diffusion fluxes, is shared by the entire density distribution along proton- 
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ospheric field lines. The density in the protonosphere (although its variation 
along field lines corresponds almost to a diffusive equilibrium) is very sensitive 
to diffusion fluxes in the transition region. 

With this understanding of the protonospheric density structure as it is 
coupled into the topside ionosphere, it is possible to outline a model of the 
plasmapause that is consistently correlated to observed ion composition in the 
topside ionosphere. 

Escape 

Figure 1 shows the field line configuration on the side of the earth [after 
Williams and Maed 1965] opposite to the sun. From this it is evident that the 
field lines become open at distances corresponding to approximately L = 6, thus 
suggesting that beyond such distances plasma escape can be effective. An esti- 
mation based on a simple escape theory will show that escape fluxes must be on 
the order of the critical flux and can therefore significantly reduce the H + con- 
centration. Thus one finds that the proton density tends to be much lower in the 
escape region than in the adjacent region where the field lines are closed and 
escape is impossible. As a result, a density gradient across field lines is built 
up, and diffusion across field lines forms a transition region between the "open" 
magnetosphere and the "closed" magnetosphere. We consider this transition 
region (shown in Figure 1 between the solid field lines) as the domain where the 
density plateau is formed that characterizes the plasmapause. 

Diffusion Across Field Lines 

From observations, in particular by Taylor et al . [1965], we can describe 
the density plateau as a rapid density decrease across field lines (from 10 3 cm -3 
at L = 4, to 10 cm -3 at L = 6 in the equatorial plane, for instance). If as assumed, 
diffusion across field lines, induced by escape, is responsible for this feature, 
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Figure 1. Field line configuration on the side 
of the earth opposite to the sun. 
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then we must consider the following physical conception. The diffusion fluxes, 
and consequently the density gradients across field lines, must be consistent 
with and coupled with the diffusion fluxes and density distributions alongj field 
lines. 

Figure 2 shows the flux pattern of a field-line section near the equator, 
also an equatorial cross-section of the density plateau. The outward flow of 
ionization across the larger field line is compensated by an inward flow across 
the shorter field line and a diffusion flux inside the tube along the field lines. 

The convergence of the field tube enhances the flux density at lower altitudes, 
and this ionization flow in the transition region between ionosphere and proton- 
osphere is responsible for the effect of the perpendicular diffusion. As pointed 
out previously, the parallel ionization flux in this region can significantly reduce 
the proton concentration in the field tube only if the flux is in the vicinity of the 
critical flux. 


If escape and diffusion across field lines are the physical processes that 
form the plasmapause and its density plateau, then the resultant ionization 
fluxes parallel to the field lines must be on the order of the critical flux. 


THEORY 

Continuity Equation 

This section presents an analytical treatment of the previously outlined 
concept. The density distribution in the protonosphere is described by the 
particle continuity equation which (under steady-state conditions) has the form: 


j r 2 3r ( r 2 Jr) + r sin 0 80 (sin 0 j^) + r 30 ]0 - 0 , 


(1) 


5 


(cm~3) 



Figure 2. Flux pattern in a field tube. Superimposed 
is the shape of the density plateau. 
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where j is the diffusion flux, and r, 0 , and 0 are spherical polar coordinates. 
For the sake of simplicity, we confine our calculation to equatorial sections of 
field tubes (as shown in Figure 2). Neglecting the 0-dependence in (1) leads then 
to 


r 2 a r ^ r2 + r 90 •*" “ 


( 2 ) 


where we introduced parallel and perpendicular diffusion fluxes 

U_l~ Jr : i„ ~ j 0 ) 

as appropriate quantities in this near-equatorial region. 

We assume that j ± is almost independent of 0 over the small region of 
concern and that j 1( = 0 at the equator, because of symmetry with respect to the 
equatorial plane. Then (2) leads to 


d_ 

dr 


< r \> = T J "e ’ 


(3) 


where here 0 is the latitude range over which we integrated (we assume 0 to 
be 30 degrees), and j is the parallel flux at this latitude (shown in Figure 2). 

As previously pointed out, the parallel diffusion flux at the bottom of the 
protonosphere can significantly affect the hydrogen-ion density. Therefore we 
transform the flux j into the corresponding flux j ( o at the bottom of the tube 
by postulating that the total flux remains constant when going down along the 
field tube. This leads to the relation 


j 


"0 


'iiq 




(4) 
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where r 0 approximately equals the earth radius. Substituting (4) in (3) gives 


_ 9 _ 

dr 


< r V = j "t> 


r 0 

e 


4 - 


3r oV /2 


( 5 ) 


V r / 

Equation 5 relates the diffusion flux across field lines to the parallel diffusion 
flux at the bottom of the protonosphere. 


Geisler and Bowhill [1965] described in some detail the dependence of 
hydrogen-ion density on the diffusion flux j in the transition region between 
ionosphere and protonosphere. At this level of our investigation we limit our 
argument by stating that the proton density in the protonosphere depends on the 
neutral-composition, oxygen-ion-density charge-exchange-rate coefficient, and, 
above all, on the transport flux j 110 . By concentrating on the relation between 
j and the proton density N 0 , we can write j |10 = j MQ (N 0 ), which, substituted in 
(5), yields 

( 6 ) 


Equation 6 describes the coupling between the perpendicular diffusion flux jj_ 
(proportional to the density gradient across field lines) and the density at the 
bottom of the field tube; this density is coupled through diffusive equilibrium 
to the proton density in the equatorial region under discussion. We can give 
(6) a more physical meaning by saying that the net diffusion flux across a field 
tube is balanced by an upward flux of ions which, in its turn, affects the density 
in the tube itself. 



Proton density is significantly affected only if the flux j,, 0 is near the 
critical flux. Then even a small variation in flux can change the proton density 
appreciably. Consequently, fluxes that can reduce the proton density so as to 
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form the described density plateau must be on the order of the critical flux j c , 
and we may therefore assume; 

j„o ~ 3c = constant . (7) 

With the assumption of (7), (6) can be integrated, which gives 



At some inner field line at a distance r, we assume that the perpendicular 
diffusion process is attenuated, and therefore the density gradient corresponds 
to diffusive equilibrium. 


Let us write for the expression 


j A =-D(r N)-^. (9) 

where D (rN) is the general expression for the altitude and density -dependent 
diffusion coefficient. This leads to 


D(rH)^-D(r,H,)^ 1 


1 iV 

J C 30 


1 



( 10 ) 


which describes the density variation across field lines in a near-equatorial 
region. 

In order to solve (10) we must know the density N 2 at some distance r 2 . 
We assume that plasma escape is what induces diffusion across field lines and 
thus forms the density plateau. The boundary value N 2 for solving (10) must 
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therefore be chosen as the density in the escape region which begins at a dis- 
tance of r 2 ~ 6r 0 . 

A very simple estimation of the escape effect leads to an approximate 
evaluation of N 2 . Following Jastrow and Rasool [1965], we find that the density 
at the escape level is less than 10 5 cm 3 . Thus H + must escape immediately 
above the charge exchange region (somewhere between 500 and 1000 km during 
minimum solar activity). With temperatures on the order of 3000 °K [Brace , 
Reddy, and Mayr , 1967] the escape flux j e is then related to the proton density 
in the ion composition transition region through 

( 11 ) 

where the subscript 2 again refers to the escape region. Obviously the upper 
limit of N 20 is determined by the highest possible flux , the critical flux j c . 
Assuming this flux to be on the order of the critical flux (assumed to be 

7 - 2-1 2—3 

j e ~ 10 cm sec ), we find a value N 20 ~ 10 cm , that is much lower than 
the diffusive equilibrium value. Conversely, if flux j e is appreciably smaller 
than the critical flux, it will not very much affect the H + density; thus N 20 will 
approach its diffusive equilibrium value, which cannot satisfy (11). Therefore, 

2 -3 

j e must be in the vicinity of the critical flux, and N 20 is on the order of 10 cm 
because of the escape of ionization. 

Assuming that the protons in the escape region are still in a state of fast 
diffusion and therefore almost in diffusive equilibrium, we find then that the 
hydrogen-ion density N 2 is on the order of N 2 ~ 10 cm 3 at r 2 = 6r 0 . After thus 
having determined the boundary conditions, (10) can be solved. In later portions 
of this paper we derive the density distributions for two kinds of diffusion mech- 
anisms and discuss the results in the light of plasmapause manifestations. 


j ~10 
J e 


X N, 


20 ! 
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Classical Diffusion 


The term "classical diffusion" refers to plasma drift across field lines, 
caused by collisions between electrons and ions. Spitzer [1956] shows that the 
perpendicular diffusion coefficient of the plasma is 

D c = 7. 2 X 10" 2 X N ■ ■ . (12) 

T n B 


where D c is the classical coefficient; N is the density; B is the magnetic field; 
and T is the plasma temperature. 

Substituting for the magnetic field in (12), the equatorial dependence 
(corresponding to the region we are concerned with) 

r 3 

and adopting from Serbu and Maier [1965] the observed temperature dependence 
as 


T 


5 X 10 2 


r 


2 



leads to 


D c = 3.5 X 10" 2 — 5 N 
r o 


Substituting this diffusion coefficient in (10) gives 


9N 2 

9r 


/jdAri 5 • 80 /Of / 3 £oV/ 2 J 3 JfcW 
\9r j[ r 5 3 c 6\rJ ^ 4 ^ 4 V 


(13) 
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which can be integrated, thus leading to 



The boundary conditions can be summarized as follows : 
For 


r 2 = 5. 5 x r 0 , N 2 = 10/cm 


(15) 


and at 



with H the density scale height of the closed, undisturbed magnetosphere. 
Applying these conditions to (14) leads to a density distribution shown in 
Figure 3 (dotted line). A closer examination of (14) reveals that the transition 
from the escape region (with a density of 10 cm -3 ) to the almost unperturbed 
"closed" magnetosphere (with its observed concentration in the order of 
10 3 cm 3 ) is extremely sharp, occurring over a distance of only a few meters. 
The diffusion coefficient used is so small that a huge density gradient is needed 
to consistently affect the proton distributions along field lines. A comparison 
of the computed density gradient with the observed density distribution by 
Taylor et al . [1966] for the equatorial region indicates a striking disagreement 
especially when we consider that the observed transition region occurs over a 
distance of almost two earth radii. 
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Turbulent Diffusion 


Failure to describe the plasmapause with the collision-type diffusion 
raises the question whether another diffusion mechanism could possibly account 
for the observed density plateau. Bohm [1943] found experimentally that the 
perpendicular diffusion became enormously enhanced for a plasma in turbulent 
state, and he verified this effect theoretically. Several authors fBraginskii . 
1965, Kauiman , 1966, and Spitzer , 1966] have since confirmed this plasma 
property. 

It is plausible that the sharp transition from the closed to the open mag- 
netosphere that is to be expected for collision-type diffusion should cause 
turbulence in the magnetosphere. We thus used Bohm's diffusion coefficient 

~ 10 3 I [Spitzer, 1956] (17) 

in our attempt to describe the plasmapause. 

If we introduce into (17) the already quoted temperature and magnetic - 
field variations, (10) leads then to the solution 



to which we apply the boundary conditions in (15) and (16). The result is shown 
in Figure 3 (solid, thick line) and reveals an acceptable agreement with the 
observed density distribution (solid, thin line). 
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It should be noted that this conformity was achieved by choosing the 
boundary level rj properly at L = 5. 5; thus, we assumed that the field lines 
become open at this distance. As shown by Taylor et al . [1966], the location 
of the plasmapause varies strongly with changing solar activity. We assume 
that this is due to variations in the field-line structure, such that there are 
variations in the location where the field lines become open. The distribution 
shown in Figure 3 should therefore be characterized by its shape and the 
approximate range of the density plateau in the transition region from the open 
to the closed magnetosphere, and not by its absolute location. 

We conclude that plasma escape in combination with turbulent diffusion 
across field lines is a valid concept for explaining the plasmapause with its 
characteristic density plateau. The resulting upward fluxes of ionization can 
evidently also account for the unexpectedly high ion-composition transition 
level at geomagnetic latitudes above 50 degrees. 

LIMITATIONS 

A rigorous treatment of the described concept would involve the simul- 
taneous solution of the nonlinear, (at least) two-dimensional and time dependent 
continuity equations for the two major ion constituents of protonosphere and 
topside ionosphere. In order to avoid this difficult task several simplifying 
assumptions and approximations were made which have to be discussed. 

The basic assumption was that the plasmapause is produced by the escape 
of ionization into the magnetic tail on the night side of the earth. Taylor's 
measurements indicate that the plasmapause similarly exists also on the day 
side where the field lines are closed; so escape cannot be directly responsible. 
Evidently, the magnetosphere is in a highly dynamic state. With proton fluxes 

7 -2 _i 

that cannot exceed the critical flux j c ~ 10 /cm sec , a protonospheric field 


15 


tube with 1 cm square basis area can only be replenished by a particle population on 
the order of 10 U during half a day. This amount is smaller by almost two 
orders of magnitude than the proton content to be expected on the basis of 
diffusive equilibrium. Thus it cannot significantly distort the distribution in- 
duced by escape; therefore, even on the day side, the typical picture of the 
plasmapause prevails. 

Furthermore, our analysis considered only the equatorial region. This 
was in order to make certain simplifying assumptions concerning the geometry 
of the field tube and the magnetic field itself, leading to an analytical solution. 

As a consequence, the ionization flux at the bottom of the protonosphere resulted 
only from a part of the field tube and is therefore somewhat low. This is so far 
not a serious limitation as other uncertainties imply deviations that could easily 
compensate for this approximation. These uncertainties concern the critical 
flux, which is not well determined because of our ignorance of the value of the 
charge-exchange coefficient; they concern Bohm's diffusion coefficient, which 
is considered as known within one order of magnitude. 

The various experimental evidences of the plasmapause are far from being 
consistent; thus the plasmapause and its form are not at all well established. 

The same is true of our concept. It is justified chiefly by the agreement between 
our theoretical result and Taylor's ion spectrometer measurements on OGO. 

We, in turn, can offer an explanation only for this type of observation. We cannot 
explain Serbu and Maier's measurements on IMP [1966], which do not suggest 
a density plateau, but rather a recovery of the plasma density in a region where 
we regard plasma escape as significant. 
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